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I The electromagnetic field in an anisotropic and inhomogeneous 



magnetodielectric is quantized by modelling the medium with two 
. independent quantum fields. Some coupling tensors coupling the elec- 

I tromagnetic field with the medium are introduced. Electric and mag- 

' netic polarizations are obtained in terms of the ladder operators of the 

medium and the coupling tensors explicitly. Using a minimal coupling 
I scheme for electric and magnetic interactions, the Maxwell equations 

and the constitutive equations of the medium are obtained. The elec- 
tric and magnetic susceptibility tensors of the medium are calculated 
in terms of the coupling tensors. Finally the efficiency of the approach 
cH I is elucidated by some examples. 
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1 Introduction 



The quantization of electromagnetic field in an absorptive dielectric, repre- 
sents one of the most and interesting problems in quantum optics, because 
it gives a rigorous test of our understanding of the interaction of light with 
matter. One of the important methods to quantize the electromagnetic field 
in the presence of an absorptive medium is known as Green function method 
p-J-[7J. In this method by adding the noise electric and magnetic polarization 
densities to classical constitutive equations of the medium, these equations 
are considered as definitions of electric and magnetic polarization operators. 
The noise polarizations are related to two independent sets of bosonic oper- 
ators. Combination of the Maxwell equations and the constitutive equations 
in frequency domain, give the electromagnetic field operators in terms of 
the noise polarizations and classical Green tensor. Suitable commutation 
relations are imposed on the bosonic operators such that the commutation 
relations between electromagnetic field operators become identical with those 
in free space. 

An interesting quantization scheme of electromagnetic field in the presence of 
an absorptive dielectric medium is based on the Hopfield model of a dielectric 
[8], where the polarization of the dielectric is represented by a damped quan- 
tum field [H]. Huttner and Barnett [TU] for a homogeneous medium and after 
Suttorp and Wubs [TT] for an inhomogeneous medium in the framework of 
the damped polarization model have presented a canonical quantization for 
the electromagnetic field inside an absorptive dielectric. This scheme is based 
on a microscopic model in which the medium is represented by a collection 
of interacting matter fields. The absorptive character of the medium is mod- 
elled through the interaction of the matter fields with a reservoir consisting 
of a continuum of the Klein-Gordon fields. In this model, eigen-operators for 
the coupled systems are calculated and the electromagnetic field operators 
have been expressed in terms of these eigen-operators. Also, the dielectric 
function is derived and it is shown that it satisfies the Kramers-Kronig rela- 
tions [12]. 

Another approach to quantizing a dissipative system is by considering the 
dissipation as a result of interaction between the system and a heat bath 
consisting of a set of harmonic oscillators [13]- [25]. In this method the whole 
system is composed of two parts, the main system and a heat bath which 
interacts with the main system and causes the dissipation of energy on it. 
In a recent approach to electromagnetic field quantization the present 
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authors have quantized the electromagnetic field in an isotropic magnetodi- 
electric [2S]. In this approach: (i) the electromagnetic field is taken as the 
main quantum system and the medium as a heat bath, (ii) The polarizability 
of the medium is defined in terms of dynamical variables of the medium, (iii) 
The polarizability and absorptivity of the medium are not independent of 
each other, as expected, this is contrary to the damped polarization model 
where polarizability and absorptivity are treated independently [IHl IH] ■ (iv) 
If the medium is both magnetizable and polarizable, one must models the 
medium with two independent collections of harmonic oscillators, where one 
collection describes electric properties and the other one describes magnetic 
properties of the medium. This scheme leads to a consistent quantization of 
the electromagnetic field in the presence of an absorptive magnetodielectric 



In the present article, the idea introduced in the previous work [2^ is gener- 
alized to the case of an anisotropic and inhomogeneous magnetodielectric. 

2 Quantum dynamics 

Electromagnetic field quantization can be achieved in an anisotropic magne- 
todielectric by modelling the medium with two independent quantum fields. 
Let us call these fields E and M quantum fields, describing the polarizability 
and magnetizability of the medium respectively. These quantum fields cou- 
ple the medium with electromagnetic field through some coupling tensors. 
The electric and magnetic polarization densities of the medium are defined 
as linear expansions in terms of the ladder operators of the E and M fields. 
The coefficients of these expansions are real valued coupling tensors. We will 
see that the electric and magnetic susceptibility tensors can be obtained in 
terms of the coupling tensors. In the following we use the Coloumb gauge 
and assume the periodic boundary conditions with no loss of generality of 
the approach. 

The electromagnetic vector potential A inside a box with volume V = 
L1L2L3 can be expanded in terms of plane waves as 



where Ufi = c\kfi\ is the frequency corresponding to the mode n, the vector 



[26]. 
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n is a triplet of integer numbers (ni, n2, ns), Y^n nicans ^ , £o is the 

ni,n2 ,n3=— oo 

permittivity of the vacuum, kfi = + + ^j^k is the wave vector 
and efix are polarization unit vectors for each n, satisfying 



enx-kn = 0. (2) 

The operators afix{t) and at^(t) are annihilation and creation operators 
of the electromagnetic field and satisfy the following equal time commutation 
rules 

[anx{t), al^y{t)] = 5n,mSxx'- (3) 

Quantization in Coloumb gauge usually needs resolution of a vector field in 
its transverse and longitudinal parts. Any vector field F{r) can be resolved in 
two components, transverse and longitudinal component which arc denoted 
by and F" respectively. The transverse part satisfy the coloumb condition 
V--F"'" = and the longitudinal component is a conservative field V x = 0. 
For a periodic boundary condition these two parts are defined as 

F^{f, t) = F{f, t)+ f (fr'V ■ F{P, t)VG{f, P), (4) 

Fll(f, t)^- I dVV • F{r', t)VG(f, P), (5) 

where 

G{r^P)^Y.T^y''-^'-''\ (6) 

n I ""n I 

is the Green function and satisfies the Poison equation 

V^G{f,P) ^ -S{f-P). (7) 

In absence of external charges the displacement field is purely transverse, 
and we can expand it in terms of the plane waves 



D{f, t) = -iso E E \/ ^ \4xit)e-''--' - anx^'^-'] enx- (8) 



n 



X=l 



2eoV 
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The commutation relations lead to the following commutation rela- 
tions between the components of the vector potential A and the displacement 
operator D 

[Ai{f,t),-D^if',t)] = ih6iif- P), (9) 

where 

5i(f - P) = i E(^^. - ^)e^'^<'-^'\ (10) 

is the transverse delta function. From (fTOj) . we see that —D plays the role 
of the momentum density of electromagnetic field. The Hamiltonian of the 
electromagnetic field inside the box is given by 



(V X Af 
2eo 2no 



A=l 



(11) 



where is the magnetic permittivity of the vacuum and we have used 
normal ordering for al^(t) and afi\(t). 

Now we include the medium in the process of quantization. For this 
purpose let the Hamiltonian corresponding to E and M quantum fields be 
denoted by He and Hm respectively. Then the medium Hamiltonian can be 
written as 



3 n+OO 

^e(t) = E E / d'qf^ddlM t)dnM t), 

3 r+oo 

Hmit) = E E / d^tuj^^M t)bnM t), (12) 

where the annihilation and creation operators dfiu{q,t), d\^{q,t), bfii,{q,t) 
and b\^{q,t) satisfy the following equal-time commutation relations 

[dnu{q,t),d^^^,{^ ,t)] = Sn,mSuu'S{q- q') , 
[bn,yiq,t),bt^,{q',t)] = 5H,rn5^u'5iq- q') ■ (13) 
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In relations fll2p uj^ is the dispersion relation of the magnetodielectric. It is 
remarkable to note that, although the medium is anisotropic in its electric 
and magnetic properties, we do not need to take the dispersion relation as a 
tensor. As discussed in |26] , we can assume a linear dispersion relation uj^ = 
c|g| with no loss of generality, but taking a linear dispersion relation simplifies 
the formulas considerably. Therefore from now on we choose the dispersion 



relation as uJq-= c\q\ where c 



is the proportionality constant. 



The basic idea in this quantization method is that the electric and mag- 
netic properties of an anisotropic magnetodielectric can be described by E 
and M quantum fields. This means that we can define the electric and mag- 
netic polarization densities of a linear but anisotropic medium as linear com- 
binations of the ladder operators of the E and M quantum fields, respectively. 
Therefore 



P,(f, = ^ E E / d'qUu;^^, r-) driu{q, t)e''^-^' + h.c 



M,{f,t) 



EE/ d^mji^^g, r) \bHu{q, t)e*^"-"' - h.c. 



n u=l 



n't 



(14) 



(15) 



where P and M are electric and magnetic polarization densities of the 
medium and 



k. 



n 



^n3 "5n3 

\'^n\ 

(16) 



In definitions of polarization densities f|T^ and f|T5l) . the real valued ten- 
sors fij{uj^, r) and gij{iOq, r), are called the coupling tensors of the electromag- 
netic field and the medium which are dependent (independent) on position r 
for inhomogeneous (homogeneous) magnetodielectrics. The coupling tensors 
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play the key role in this method and are a measure for the strength of the 
polarizability and magnetizability of the medium macroscopically. We will 
see that the imaginary parts of the electric and magnetic susceptibilty in 
frequency domain can be obtained in terms of these coupling tensors. Also, 
explicit forms for the noise polarization densities can be obtained in terms 
of the coupling tensors and the ladder operators of the medium. The cou- 
pling tensors are common factors in the noise densities and the electric and 
magnetic susceptibilities, and so the strength of the noise densities are de- 
pendent on the strength of the electric and magnetic susceptibility. It can 
be shown that for a non absorptive medium, the noise densities tend to zero 
as expected and this quantization scheme reduces to the usual quantization 
in such media. 

A consistent quantization scheme must lead to the correct equations of 
motion of the system and the medium. These equations are macroscopic 
Maxwell and constitutive equations of the medium and we will see that these 
equations can be obtained from the Heisenberg equations using the total 
Hamiltonian defined by 



d r 



[D(f, t) - P(f, t)f ^ (v.Anr.t) _ ^ ^ ^^^^ . ^1 



2Sn 



(17) 



2.1 Maxwell equations 

Using the commutation relations ^ the Heisenberg equations for the vector 
potential A and the displacement field D are 

(18) 

^ . hH,Dm\ - ^ii^iiiM - V X M(f,t), (19) 
ot n /io 

where is the transverse component of P. The transverse electrical field 
i?-*-, magnetic induction B and magnetic field H are defined by 

dA B 

E^ = - — , B = VxA, H= M. (20) 

dt Ho 
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Using these recent relations, fllSp and fll9l) can be rewritten as 

D = soE^ + P^, (21) 
dD 

- = Vx«, (22) 

which are the definitions of the displacement field and the macroscopic Maxwell 
equation, in the presence of a magnetodielectric, respectively. 



2.2 Constitutive equations of the medium 

A magnetodielectric subjected to electromagnetic field can be polarized and 
magnetized in consequence of interaction of the medium with the field. The 
macroscopic electric and magnetic polarizations is related to electric and 
magnetic fields, respectively by the constitutive equations of the medium. 
Therefore a quantization scheme must be able to give the constitutive equa- 
tions in the Heisenberg picture. In this section by applying the Heisenberg 
equations to the ladder operators of the medium we find the correct consti- 
tutive equations of the medium. 

the The time evolution of the operators dnu{q,t) and bfii,{q,t) can be 
obtained from the commutation relations (fT3l) and the Hamiltonian (fTTll as 
follows 

t ~ 

dnAq,t) = -[H,dri^{q,t)] = 

-luj^nM + [ d'r'e-'^--^^'Uuj^, r')E\r' , t)vl, (23) 

h 

-iLo^nM t) + -^ [ d^?e-''^--^\,,{uj^, r')B\?, t)si^. (24) 
nyV Jv 

It is easy to show that these equations have the following formal solutions 



dHMt)=dnM^y''' + 



dt'e'^"f(*-*') / d'r'e 



hVv Jo 



(25) 
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hVV Jo 



V 



Substituting (!25l) in (fT4l) and ( !26l) in ( fTSll we obtain the macroscopic con- 
stitutive equations of the anisotropic polarizable and magnetizable medium, 



1*1 



F(r, t) = P^ir, t)+eo dt'x\r, \t\ - t')E{r, ±t'), 

Jo 

1 /-I*! 

M(r, t) = MN{r, t) + — dt'x"^ir, \t\ - t')B{f, ±t') 

/in Jo 
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1 /-I*! 

where the upper (lower) sign corresponds to t > (t < 0) and E 
is the total electrical field. 
The memory tensors 



(27) 
(28) 



dA _ Pi 

dt eo 



X"(r,t) 



f j^J,^duu\fn{u,r)smu;t 


^/o°°da;a;2(^^t)(a;,f)sinu;t 







t > 0, 

t < 0, 
t > 0, 

t < 0, 



(29) 



(30) 



are called the electric and magnetic susceptibility tensors of the mag- 
netodielectric, respectively and /*, denote the transpose of the coupling 
tensors /, g. If we are given a definite pair of tensors x'^(^i t)y X™'(^) t) which 
are zero for t < 0, then we can inverse (129|) and (l30l) and obtain the corre- 
sponding tensors (//*) and (gg^) as, 



(//*)(^,r-) 



C t^^X^r, t) sincut = ^Im \r{ 



[r,LO, 



u > 0, 
w = 0, 
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(31) 



{gg'){uj,r) = 



Tlc? too 



S^dtr^ir^t) smut 



■Im 



iu > 0, 
uj = 0, 



(32) 

where x'^(r, a;) and x™'(r, a;) are the susceptibihty tensors in the frequency 
domain. The operators Pjy and Mn in and are the noise electric 
and magnetic polarization densities 



Pm{r,t) 



3 

EE 

n 1^=1 



h.c. 



(33) 
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si 



(34) 



These noises are necessary for a consistent quantization of the electro- 
magnetic field in the presence of an absorptive medium. 

From fl3T|) and fl32l) it is clear that for a given pair of the susceptibility 
tensors and there are infinite number of coupling tensors / and g 
satisfying the equations fl3ip and ([32D- In fact for a given pair of and 
X"* if the tensors / and g satisfy equations (I3T1) and (132|) . then the coupling 
tensors fA and gA, for any orthogonal matrix {AA^ = 1), are also a solution. 
Certainly this affect the space-time dependence of the noise polarizations and 
therefore the space-time dependence of the electromagnetic field operators, 
but all of these choices are equivalent. This means that the various choices of 
the coupling tensors / and g satisfying fl3Tl) and fl32l) . for a given pair of the 
susceptibilities and x™, do not affect the commutation relations between 
the field operators and hence the physical observables. This becomes more 
clear if we compute the commutation relations between the components of 
the Fourier transform of the noise polarizations 



10 



heo 



X?.(r,a;) S{f - r')6{uj ~ u'), 

L — J 

—Im x'^(r, u;) S{r' — r')6{uj — uj'). 



71 



h 



(35) 



These relations are generalization of those in reference [27]. For a given 
pair of and x™, various choices of the coupling tensors / and g satis- 
fying the relations fl3T|) and fl32|) . do not affect these commutation relations 
and accordingly the commutation relations between the electromagnetic field 
operators. Hence, all of the field operators which are obtained by using a 
definite pair of the susceptibilities and x™, with different coupling tensors, 
satisfying fl^ and (1521) . are equivalent. 

3 Solution of Heisenberg equations 

The Maxwell and constitutive equations of the medium constitute a set of 
coupled equations. In this section we solve them in terms of their initial con- 
ditions using the Laplace transformation technique. For any time-dependent 
operator g{t) the forward and backward Laplace transformation of g{t) are 
defined by 



respectively. Carrying out the forward and backward Laplace transfor- 
mation of the Maxwell equation fl22l) and the constitutive equations (I2T]) . 
(EZD and (EH]) we find 



where e{f, s) = 1 + x'^if, s) and x"^{f, s) are the Laplace transformations 
of the electric permeability tensor and magnetic susceptibility tensor of the 
medium, respectively and 




(36) 



V X V X E^^\f, s) + fioeoS^e{f, s)&^\r, s) 

V X x'"(r, s)V X &^\r, s) = jf^\r, s), 



(37) 
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{r, s) = ±V X S(f, 0) - fios^Pi^\f, s) T 
/iosV X Mlj\f, s) T V X x"^(r, s)B{r, 0) + fiosD{f, 0), (38) 



is the forward and backward Laplace transformation of the noise cur- 
rent where upper(lower) sign corresponds to J-^(r, s) {J-^{f,s)). The wave 
equation fl37j) can be solved using the Green tensor method |27j. To see the 
space-time dependence of electric field more explicitly, let us consider a ho- 
mogeneous but anisotropic bulk medium. In this case by expanding E-^'^{f, s) 
and Jf'^{f, s) in terms of plane waves as 



(39) 



and inserting this expansions in the wave equation fl37p we obtain 
A(^n, s)E^'^{kn, s) = f^^{K-t, s), 

Aij{kn, s) = ~eif,u£af3j [Sua - xTai^)] k'tk^ + ^QeQS^eij{s), (40) 

where Si^^ is the Levi-Civita symbol, we can use the expansions ([I]), ([8]), (133!) 

and (IMI) to obtain the operator J ' (fc^j, s) in terms of the ladder operators of 
the electromagnetic field and the magnetodielectric medium. Finally using 
(15^ and (1401) after some elaborated calculations we obtain the space-time 
dependence of the electric field in terms of the ladder operators of the field 
and medium as 



2cV 



n A=l 
3 

^° EE [d'q[^t,{uJ,^kH,t)dHAme''^-''+h.c. 
±^EE / d'<l [CfM''^H,t)bHMO)e'''-'+h.c: 



V 
fJ-O 



(41) 
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where the upper (lower) sign corresponds to t > (t < 0) and ?7j^(A;^, +t), 

Vij{kn,-t), itji.^q^kfi,+t), iij{uj^,kH,-t), Ctji.^qM^+t) and Ciji.^qM^-i) 
for t > are given by 



S It ILug 



L 1 <^ Au\kri, s)— \ eia(3k^gpj{uj^), 




(42) 



and L^^{h{s)} denotes the inverse Laplace transformation of h{s) and 
is the inverse of the matrix A. 
Example 1: 

In the first example we show that in the absence of any medium this quanti- 
zation scheme reduces to the usual quantization in the vacuum. In free space 
the electric and magnetic susceptibility tensors are zero and from flHT]) and 
(!32|) we deduce that the coupling tensors / and g are also zero. Therefore 
from (go]), gl]) and (HID one finds 



which is the electric field in the free space. So in this case, quantization 
of electromagnetic field is reduced to the usual quantization in the vacuum 
as expected. 
Example 2 : 

Take the susceptibility tensors x"^ and as follows 




(43) 




A 



< t < A, 







otherwise. 




otherwise 



< t < A, 
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where Xo(r) and Xo'i^ are some time independent but position dependent 
tensors and A is a real positive constant, using (1ST]) and (152]) we find 



{gg'){uj,r) 



(44) 



and from (ETj) and (EH 



P(f, t) = P^{f, t) + x^(r-)^ /'*' dt'E{f, ±t'), 

M(r, t) = M^ir, t) + Xoi^^ ^ (4^) 

where PN{f,t),M]\[{f,t) are the noise polarization densities correspond 
to a pair coupling tensors / and g satisfying (1^4]) . In the limit A ^ 0, from 
(jH]) we deduce that the coupling tensors and therefore the noise polarization 
densities defined by ([23]) and (IM]) tend to zero. In this case the constitutive 
equations (115]) are 



(46) 

and the wave equation (157j) becomes 



V X V X P^\f, s) + fioEos^ [1 + Xoir)] Ef^\r, s) - 

V X x^(f)V X E^'\f, s) = ,P'\f, s), (47) 

where the noise current density (138]) is 
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P'\r, s) = ±V X B{r, 0) ^ V x Xo{r)B{f, 0) + fiosD{f, 0). 

(48) 

We see that the noise operators have vanished. This is because in the 
hmit A — > 0, the absorption coefficients tend to zero. The solution of the 
wave equation (H7j) can be expressed in terms of the Green tensor as 

Ef'^{f, s)= J d^r'G{f, r', s)jf'\r', s), (49) 
where the Green tensor satisfies the equation 

V X V X G{r, P, s) + /io^os^ [1 + Xoir)] G{r, ?, s) - 

V X x[r(^ V X G{f, ?, s) = S{f- P), (50) 

together with some boundary conditions. These boundary conditions 
guarantee the continuity of the tangential component of electric field and the 
normal component of magnetic field at some surfaces where the susceptibil- 
ities of the medium become discontinuous. For an anisotropic homogeneous 
medium using PUI) . fHTl) and fH2|) we can write the electric field as follows 



n A=l 



2cV 



7]tJkH,t)anx{0)e'''"' + h.c. 



(51) 



where 



r)fj{kH, ±t) = L ^ |Ai/(A;^, s) {is ± loh) Sij ± ujHeifiu£af3jk'ik§{x'S')^a } , 
Au{kn, s) = -Si^^Eo^pi [5uc, - (Xo%«] + /"o^os^ (1 + {xl)ii) ■ (52) 

This example shows that the present quantization scheme is reduced to 
the usual quantization in a nonabsorptive medium. 
Example 4: A simple model for electric susceptibility tensor 

If we neglect the difference between local and macroscopic electric field for 
substances with a low density, then the classical equation of a bound atomic 
electron in an external electric field for small oscillation can be written as 
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Xi + 270;, 



E^{t), 



m 



1,2,3, 



(53) 



where E{t) is the electric field in the place of the atom and the magnetic 
force has been neglected in comparison with the electric one and 7 is a 
damping coefficient. We have assumed that for sufficiently small oscillations 
the ith component of the force exerted on the bound electron by nucleus, 
can be expressed as a linear combination of the coordinates of the electron 
with constant coefficients Kij. Therefore in this simple model the motion of 
the bound electron is as a forced anisotropic harmonic oscillator. Let E_{^) 
and f{u) be Fourier transforms of the electrical field E{t) and position f{t) 
respectively. From ( l53l) we find 



r(uJi 



e 

m 



V + 2z7u;)l + K 



E(uj) 



(54) 



Now let there be molecules per unit volume of the medium with z 
electrons per molecule. We assume that the damping coefficient (7) and the 
tensor K are identical for each electron. Then for the Fourier transform of 
the electric polarization density we find 



P(r 



,UJ 



m 



-UJ 



+ 2z7Cj)l + K K{r,u)) 



(55) 



From fl55|) we find the electric susceptibility tensor of the medium in 
frequency domain 



From fl3Tl) we have 



{-u? + 2z7u;)l + K 



(56) 



(//*)( 



§^270; + Kf + 472u;2ll " a; 7^ 0, 



(57) 











In the special case 7 = 0, the anisotropic dielectric substance is a nonab- 
sorptive one and this relation for oj 7^ becomes 



(//*)0 



ViXi yf ViZi 



(58) 
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where Ui, are eigenvalues of the tensor K corresponding to eigenvectors 

Ri = \ Hi \i = 1, 2, 3). In this case the couphng tensor / is nonzero only 

for frequencies c^i, ^^2, i^s,- These frequencies are the resonance frequencies 
of the equation (!53|) . Therefore when 7 = 0, that is for a non absorptive 
medium, the coupling tensor / and therefore the noise electrical polarization 
density is equal to zero except in resonance frequencies, where the energy of 
electromagnetic field is absorbed by the oscillator. This example explicitly 
shows that this quantization scheme is also applicable to anisotropic disper- 
sive but non absorptive media. 



4 Summary 

The electromagnetic field quantization in the presence of an anisotropic mag- 
netodielectric is investigated consistently by modelling the magnetodielectric 
with two independent quantum fields namely E and M quantum fields. For 
a given pair of the electric and magnetic susceptibility tensors ^'^^ X^i 
we have found the corresponding coupling tensors / and which couple 
electromagnetic field to E and M quantum fields respectively. The explicit 
space-time dependence of the noise polarizations are obtained in terms of the 
ladder operators of the medium and the coupling tensors as a consequence 
of Heisenberg equations. Maxwell and constitutive equations are obtained 
directly from Heisenberg equations. In the limiting case, i.e., when there is 
no medium, the approach tends to the usual method of quantization of the 
electromagnetic field in vacuum. Also when the medium is a non absorptive 
one, the noise polarizations tend to zero and in this case the approach is 
equivalent to the previous methods, as expected. 
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